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Abstract New metastable massive particles with elec- 
tric and colour charge are features of many theories 
beyond the Standard Model. A search is performed for 
long-lived gluino-based i?-hadrons with the ATLAS de- 
tector at the LHC using a data sample corresponding to 
an integrated luminosity of 31 pb _1 . We search for evi- 
dence of particles that have come to rest in the ATLAS 
detector and decay at some later time during the peri- 
ods in the LHC bunch structure without proton-proton 
collisions. No significant deviations from the expected 
backgrounds are observed, and a cross-section limit is 
set. It can be interpreted as excluding gluino-based R- 
hadrons with masses less than 341 GeV at the 95% C.L., 
for lifetimes from 10~ 5 to 10 3 seconds and a neutralino 
mass of 100 GeV. 



1 Introduction 

The search for exotic massive long-lived particles (LLPs) 
is an important component of the early data exploration 
program of the Large Hadron Collider (LHC) experi- 
ments. LLPs are predicted in supersymmetry (SUSY) 
models, such as split-SUSY [TJ[2] and gauge-mediated 
SUSY breaking [3 , as well as other exotic scenarios, e.g. 
universal extra dimensions [3]. If LLPs carried colour 
charge, as predicted in a variety of Standard Model 
(SM) extensions [S], they would hadronise with quarks 
and gluons to form colour-singlet states before interact- 
ing with the detector. i?-parity conserving split-SUSY 
presents such a model since the gluino (g) decay — g — > 
gXi, or g — > qq~Xi — is suppressed by small coupling to 
the gravitino LSP or very large squark (q) mass 6 . The 
possibility of direct pair production of coloured parti- 
cles through the strong nuclear force implies a large pro- 
duction cross-section at the LHC. The search described 



herein is particularly sensitive to such a scenario, and 
results are interpreted in the context of split-SUSY. 

During the hadronisation process, long-lived gluinos 
bind with SM quarks and gluons from the vacuum to 
produce "i?-hadrons"Q — a new heavy composite state [7] . 
When i?-hadrons scatter via the nuclear force, they can 
transform their internal state by exchanging partons 
with the detector material (e.g. gqq —> gqqq). Depend- 
ing on the mass hierarchy of i?-hadron states and the 
matter interaction model, these particles traverse the 
detector flipping between electrically neutral, singly- 
and doubly-charged states. Several properties of the R- 
hadron, such as mass and lifetimqj, are primarily de- 
termined by properties of the constituent gluino. Other 
properties, such as hadronic matter interactions, are 
dominated by the constituent quarks and gluons. 

Since i?-hadrons may be produced near threshold 
in LHC collisions, they are expected to be slow (/3 sig- 
nificantly below 1) and have large ionisation energy 
loss. At LHC energies, some fraction (typically several 
per cent) of i?-hadrons, produced with low kinetic en- 
ergy, would lose sufficient energy to come to rest inside 
the dense detector materials of the ATLAS calorime- 
ter. These stopped i?-hadrons may have lifetimes span- 
ning many orders of magnitude, and may decay with 
significant delay after the collision that created them. 
It is possible to more easily detect these decays by 
searching for calorimeter energy deposits during the 
so-called empty bunch crossings [8] when there is less 
background, as described in detail in Section [3] Since 
this sample has a very low proton-proton collision rate, 
the dominant backgrounds to this search are cosmic 



1 The term i?-hadron has its origin in the _R-parity quantum 
number in supersymmetry theories. 

Lifetime refers to the gluino decays, not the _R-hadron tran- 
sition from one hadronic state to another. 
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ray muons, beam-related backgrounds and instrumen- 
tal noise. Similar searches for out-of-time decays have 
previously been performed by other experiments [9lll0j . 

This analysis complements previous ATLAS searches 
for long-lived particles [TTl[T2j which are less sensitive 
to particles with initial ft -C 1. By relying primarily 
on calorimetric measurements, this analysis is also sen- 
sitive to events where i?-hadron charge-flipping may 
make reconstruction in the inner tracker and the muon 
system impossible. A potential detection of stopped 
i?-hadrons could also lead to a measurement of their 
lifetime and decay properties. Moreover, the search is 
sensitive to any potential new physics scenario produc- 
ing large out-of-time energy deposits in the calorime- 
ter with minimal additional detector activity. The data 
analysed in this Letter were recorded by the ATLAS 
experiment between April and November 2010, exploit- 
ing proton-proton collisions at a centre-of-mass energy 
of 7TeV. 



2 The ATLAS detector 

The ATLAS detector 02] consists of an inner track- 
ing system (ID) surrounded by a thin superconduct- 
ing solenoid, electromagnetic and hadronic calorime- 
ters and a muon spectrometer (MS). The ID consists 
of pixel and silicon microstrip detectors, surrounded by 
a transition radiation tracker. The calorimeter system is 
based on two active media for the electromagnetic and 
hadronic calorimeters: liquid argon (LAr) in the inner 
and forward regions, and scintillator-tiles (TileCal) in 
the outer barrel region. The absorber materials are lead 
and either steel, copper, or tungsten. The calorimeters 
are segmented in cells which have typical size 0.1 by 0.1 
in r]-(f) spac^l in the TileCal section. The MS, capable 
of reconstructing tracks within \rj\ < 2.7, uses toroidal 
bending fields generated by three large superconducting 
magnet systems. 

Jets are constructed using the anti-fct jet algorithm [2] 
with a radius parameter (in rj-(j) space) set to R = 0.4, 
which assumes the energetic particles originated from 
nominal interaction point. This assumption, while in- 
correct, still accurately quantifies the energy released 
from the stopped i?-hadron decays occurring in the 
calorimeter; the associated systematic uncertainty is 

3 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The z-axis 
points from the IP to the centre of the LHC ring, and the 
j/-axis points upward. Cylindrical coordinates (r, (f>) are used 
in the transverse plane, <j> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as r\ = — lntan(#/2). 



discussed in Section [7] The inputs to the jet algorithm 
described in this paper are calorimeter energy deposits. 
Jet energy is quoted without applying a hadronic cal- 
ibration; only jets with energy above 10 GeV are con- 
sidered in this search. ATLAS jet reconstruction algo- 
rithms are described in more detail elsewhere |15j . 

3 LHC bunch structure and trigger strategy 

The LHC accelerates two counter-rotating proton beams, 
each divided into 3564 25 ns bunch slots. When protons 
are injected into the LHC not every slot is filled. In 
late 2010 running, slots that were filled typically had 
10 11 protons. Unfilled slots could contain protons due 
to diffusion from filled slots, but this was typically be- 
low 5 x 10 8 protons per slot. 

The filled and unfilled slots can be combined to 
make three different "bunch crossing" scenarios. A 'paired 
crossing consists of a filled bunch from each beam col- 
liding in ATLAS and is by far the most likely to produce 
i?-hadrons. An unpaired crossing has a filled bunch from 
one beam and an unfilled slot from the other. Finally, 
in an empty crossing the slots from both beams are un- 
filled. Empty crossings typically had a proton-proton 
collision rate less than one-millionth the rate in paired 
crossings [15]. In the search sample, described in Sec- 
tion [SJ there are approximately 350 paired bunch cross- 
ings organized in groups of eight. Inside a group, paired 
bunch crossings occurred every 150 ns or longer, while 
the gaps between groups are longer. Most ATLAS anal- 
yses use data collected from the paired crossings; this 
analysis instead searches for physics signatures of LLPs 
in the empty crossings. This is accomplished with a set 
of dedicated low-threshold calorimeter triggers which 
may fire only in the empty or unpaired crossings where 
the background to this search is much lower. 

ATLAS has a three-level trigger system consisting of 
one hardware and two software levels [T7]. Signal can- 
didates for this analysis are collected using a hardware 
trigger requiring localised calorimeter activity with a 
10 GeV transverse energy threshold. This trigger could 
fire only during an empty crossing at least 125 ns af- 
ter the most recent paired bunch crossing. By waiting 
to collect data only during the empty crossings many 
signal decays are lost; however, this provides a sample 
nearly free of all collision backgrounds. Although the 
i?-hadrons decay at a time randomly distributed with 
respect to the bunch-crossing clock (which has a 25 ns 
period), this does not cause any loss of efficiency as the 
calorimeter response is longer than 25 ns. A sideband 
region to study beam-halo muons is collected with a 
similar trigger that fired in the unpaired crossings. Both 
samples are collected with only a hardware-level trigger 



3 



and without further requirement at the higher trigger 
levels. 



4 Simulation of i?-hadrons 

Monte Carlo simulations are used primarily to deter- 
mine the reconstruction efficiency and stopping fraction 
of the i?-hadrons, and to study associated systematic 
uncertainties on the quantities used in the selections. 
The simulated samples have gluino masses in the range 
200-600 GeV, to which the present analysis is sensitive. 
Furthermore, a detailed simulation incorporating gluino 
generation, stopping and decay steps provides a mod- 
ular approach for signal production. The Pythia pro- 
gram 1 18] is used to simulate gluino-gluino pair produc- 
tion events. The string hadronisation model (TH], incor- 
porating specialised hadronisation routines [5] for R- 
hadrons, is used inside Pythia to produce final states 
containing two i?-hadrons. 

Table 1 Signal Monte Carlo stopping fractions for various 
gluino mass values with their statistical uncertainties. Each 
efficiency is the number of _R-hadrons that stop anywhere in 
the detector divided by the number of produced fl-hadrons 
events. Some B-hadrons stop in parts of the detector where 
there is little sensitivity to detect their decay. This is ac- 
counted for in the reconstruction efficiency in Table [2] 



mg (GeV) | Generic (%) \ Regge (%) 



200 


11.5±0.7 


5.7±0.5 


300 


13.1±0.7 


4.7±0.5 


400 


12.4±0.7 


6.5±0.5 


500 


13.2±0.7 


4.5±0.5 


600 


11.8±0.7 


5.6±0.5 



To compensate for the fact that i?-hadron scatter- 
ing is not strongly constrained by SM analogues, the 
simulation of i?-hadron interactions in matter is han- 
dled by a special detector response simulation [50] using 
Geant4 [2"Tll2"2"] routines based on two rather different 
scattering models with different sets of assumptions: the 
Generic [2"Ull2"5] and Regge [2"4"ll2"5] models. Each model 
makes different assumptions about the i?-hadron nu- 
clear cross-section and mass spectra of various internal 
states. Briefly, the phenomenologies of the two models 
are described as follows: 

Generic : Limited constraints on allowed stable states 
permit the occurrence of doubly charged i?-hadrons 
and a wide variety of charge- exchange scenarios. 
The scattering model is purely phase space driven. 
This model is chosen as the nominal signal model 
for gluino i?-hadrons. 



Regge : Only one (electrically neutral) baryonic state 
is allowed. The scattering model employs a triple- 
Regge formalism. This is expected to result in a 
lower stopping fraction than the Generic scenario. 

If an i?-hadron comes to rest in ATLAS, its location is 
recorded. TableQ]shows the probability for an i?-hadron 
to stop as a function of the generated gluino mass for 
the two models considered. The stopping fraction shows 
no significant dependence on the gluino mass within the 
available simulation statistics. 

These stopping locations are used as input into a 
second step of Pythia where the decays of the R- 
hadrons are simulated. Different models allow the gluinos 
to decay via the radiative process, g — > gx\, or via 
g — > qqxl- The reconstruction efficiencies are compa- 
rable for both decay modes however, the results are 
interpreted assuming a 100% branching ratio to gxl 
in accordance with previous results [91110] . In all simu- 
lations the neutralino mass, m^o , is fixed to 100 GeV. 
Interactions of the decay products with the detector are 
simulated with Geant4. These events then follow the 
standard ATLAS reconstruction which outputs signal 
candidates. 
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Gluino Lifetime (seconds) 

Fig. 1 Percentage timing efficiency as a function of gluino 
lifetime. The region 10 -5 to 10 3 seconds has the highest ef- 
ficiency, approximately 37%. For shorter lifetimes many R- 
hadrons decay in paired bunch crossings while for longer life- 
times many decay when ATLAS is not taking data. The solid 
line is calculated with the assumption that no fl-hadron sur- 
vived from one run to the next (allowing a per run analysis for 
shorter lifetimes). The dashed line is calculated after averag- 
ing over bunch structure but does allow fl-hadrons to decay 
in a run they are not produced in. 

An additional inefficiency arises because only a frac- 
tion of the stopped i?-hadrons decay in an empty cross- 
ing, while ATLAS is taking data. This inefficiency is a 
function of the gluino lifetime, and is calculated from 
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a detailed knowledge of the bunch structure across dif- 
ferent LHC fills (relevant for short lifetimes) and the 
ATLAS data acquisition schedule (relevant for long life- 
times) . The timing efficiency versus lifetime is shown in 
Figure [TJ the plateau corresponds to 37% efficiency for 
lifetimes from 10~ 5 to 10 3 seconds. These lifetimes are 
long enough to survive until an isolated empty cross- 
ing which may be many microseconds after the produc- 
tion event. However, these lifetimes are short enough 
so that most of the i?-hadrons do not decay outside of 
an ATLAS data-taking run, which is typically ten to 
twenty hours. 



6 Candidate selection 

In addition to the signal candidates, several background 
processes survive the trigger: beam-halo muons, proton- 
proton collisions, proton-gas collisions, calorimeter noise 
and cosmic ray muons. The first three of these processes 
arise from the non-zero proton population in the empty 
bunches, as described in Section [3] Beam-halo muons 
are produced when a stray 3.5 TeV proton collides with 
a collimator or residual gas many metres away from 
ATLAS. The subsequent shower decays or is blocked 
by shielding, but occasionally a muon survives to tra- 
verse the detector. 

The search and sideband samples are divided into 
two exclusive channels based on the number of recon- 
structed jets with energy over 10 GeV, Nj cts - Lower en- 
As discussed in Section^ data are collected with calorimetef^gy ^-hadron decays tend to be reconstructed as a sin- 
based triggers in the empty and unpaired crossings. To S le J et ' but ver y energetic ones may be reconstructed as 



5 Data samples 



construct two sidebands and a signal search sample, the 
data from empty bunch crossings are separated into 
early, middle and late 2010 subsets, respectively. The 
total integrated luminosity in paired crossings concur- 
rent with these three empty-crossing samples is 34 pb -1 
[IB] . Each sample has a different signal to background 
ratio due to the rapidly increasing instantaneous lumi- 
nosity of the LHC during 2010. 

To predict the number of cosmic ray muons — the 
dominant background — both the amount of time and 
varying number of empty bunches per beam revolution 
in each data-taking run must be accounted for. The 
trigger live-time (in units of crossing-hours) is quanti- 
fied by taking the product of the number of hours the 
trigger is active and the number of empty bunch cross- 
ings per beam revolution in that LHC filling scheme. 

The first data sample, the background sample, had 
low beam current and corresponds to a trigger live-time 
of 1.07xl0 6 crossing-hours and 0.26 pb -1 of integrated 
luminosity; it is used to estimate the background from 
cosmic ray muons. A second sample, the control sam- 
ple, corresponds to 0.55 x 10 6 crossing-hours, and an in- 
tegrated luminosity of 2.6 pb -1 . This sample has dif- 
ferent beam conditions, and acts as a cross-check for 
the cosmic ray muon background estimate. The signal 
search is performed in the search sample, which corre- 
sponds to 0.32xl0 6 crossing-hours. Data collected in 
the search sample reflect an integrated luminosity of 
31 pb -1 with a peak instantaneous luminosity of ap- 
proximately 2.1 x 10 32 cm~ 2 s _1 . The probability of pro- 
ducing an i?-hadron in either the background or control 
sample is small compared to the search sample. 

Finally a small sample from the unpaired crossings 
is used to quantify the background induced by beam- 
halo muons. 



multiple jets due to their larger spatial extent. The two 
channels, denoted single jet and multi-jet, differ primar- 
ily in their respective noise rejection criteria, labelled as 
"jet cleaning" . Without requirements on a second jet 
to reject background, the single jet channel, Nj ets = 1, 
uses a more demanding calorimeter selection. Selection 
criteria described below that are not specifically men- 
tioned in the event yield summary (Table [2]) fall under 
the jet cleaning label for either the single or multi-jet 
channel. 



6.1 Criteria common to both channels 

All events in the three empty-crossing data samples, de- 
scribed in Section[5l must pass basic beam, detector and 
data quality requirements. Events with proton-proton 
collisions or proton-gas collisions, which happen at a 
very low rate in empty crossings, are rejected by re- 
quiring no reconstructed tracks in the ID. Additionally, 
events where the missing transverse momentum is less 
than half of the leading jet transverse momentum are 
vetoed. The criterion on missing transverse momentum 
removes the rare collision events since they tend to pro- 
duce jets with balanced transverse momentum. A signal 
event, however, produces only one localised energy de- 
posit, since even if both i?-hadrons stop in the detector 
they are likely to decay at times separated by more 
than 25 ns. The trigger system would identify these de- 
cays as separate events because its timing resolution is 
significantly smaller than 25 ns. 

To reduce the contribution from calorimeter noise, 
the most energetic jet must have an energy greater than 
50 GeV, and be in the central region (\rj\ < 2.2). Fur- 
thermore, ngo, the minimum number of calorimeter cells 
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Table 2 Event yield at various levels of the selection for the single jet (top) and multi-jet (bottom) channels. The yields of 
the background and control samples are scaled by detector live-time to the search sample. The background and control samples 
agree well in the bottom three rows for both channels, showing the live-time scaling works after detector noise is removed (via 
jet cleaning). Agreement is not expected before jet cleaning due to the sporadic nature of calorimeter noise across 2010. Since 
the rate of beam-halo muons scales with beam current it only contributes significantly to the search sample as discussed in 
Section [8] Only statistical uncertainties are shown in the table. 



Selection criterion I Background sample I Control sample I Search sample 



Data Quality 


1119200 ± 600 


1038600 ± 800 


1579949 


|Jet ri\ < 1.2 


875600 ± 500 


739600 ± 700 


1094196 


N jeta = 1 


71800 ± 200 


72800 ± 200 


1089374 


Jet Cleaning 


5860 ± 40 


5860 ± 60 


5615 


Muon Segment Veto 


4 ± 1 


5 ± 2 


9 


Jet Energy > 100 GeV 


0.3 ± 0.3 


0.6 ± 0.6 





Data Quality 


1119200 ± 600 


1038600 ± 800 


1579949 


|Jet rj\ < 2.2 


1093300 ± 600 


1012400 ± 800 


1507017 


K iVjots <10 


22480 ± 80 


22700 ± 100 


24902 


Jet Cleaning 


8650 ± 50 


8310 ± 70 


8036 


Muon Segment Veto 


1 ± 0.6 


4 ± 2 


3 


Jet Energy > 100 GeV 


0.6 ± 0.4 


0.6 ± 0.6 


1 



that collectively contain 90% of the jet energy, must be 
greater than three. 

Cosmic ray and beam-halo muons are vetoed by re- 
quiring the absence of a muon segment anywhere in 
the MS. A muon segment is formed when there is cor- 
related activity in a single detector station of the MS, 
and serves as input to more sophisticated track find- 
ing algorithms. This presents the common selection in 
the different jet channels and is used to study the data 
sample. In the final search only events with leading jet 
energy over 100 GeV are considered. This threshold op- 
timises the search reach while maintaining sensitivity 
to lighter gluino models. 



6.2 Single jet criteria 

To remove effects from the noisier calorimeter endcaps, 
the leading jet must be well within the barrel region 
(\r]\ < 1.2). Jets must have at least half of their energy 
in the TileCal (/tiio > 50%) to reject beam-halo muons, 
which occur more frequently closer to the beam pipe. 
The energy weighted transverse size of the jet, in r\ — (\> 
space, must be greater than 0.04. 



6.3 Multi-jet criteria 

The requirement of a second jet in the multi-jet chan- 
nel, 1 < A/jets < 10, strongly suppresses the noise and 
beam-halo contributions, allowing for looser jet clean- 
ing cuts. The energy of the sub-leading jet must be at 
least 15 GeV to ensure it is well reconstructed; requiring 
Aj c t s < 10 vetoes several events with bursts of calorime- 
ter noise. In this channel, there is no jet transverse size 



requirement; however, one of the two leading jets must 
have /riic > 10%. 



6.4 Data sample consistency 

A very small fraction of the cosmic ray muons incident 
on the ATLAS detector deposit sufficient energy in the 
calorimeters to mimic the i?-hadron signal. This is the 
dominant background in all data samples prior to ap- 



plying the muon segment veto. Figures 2(a) and 2(b) 
show the jet energy distributions, scaled to live-time, of 
the events selected in the single and multi-jet channels, 
respectively. The distributions show the event popula- 
tions after applying the jet cleaning criteria but prior 
to the muon segment veto. Agreement among the three 
empty-crossing data samples is found, demonstrating 
that cosmic ray muons dominate the background at this 
point. 



6.5 Signal efficiency 

The same selection criteria are applied to both the data 
and the signal sample; their affect on signal efficiency is 
calculated using the simulation. There is an additional 
small signal loss when a muon segment is reconstructed 
in the same event as the signal decay. Muon segments 
due to cosmic rays or noise would cause the event to 
be vetoed but they are not included in the simulation, 
unlike those due to signal decays. Using data from the 
empty bunches acquired with a random trigger, it is es- 
timated that 93% of decays would pass the muon seg- 
ment veto and remain in the search sample. 
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Fig. 2 Distributions of jet energy for the single jet (a) and 
multi-jet (b) channel. Events are plotted for the three empty- 
crossing data samples considered in the analysis, scaled to 
the detector live-time of the search sample, prior to applying 
the muon segment veto. Cosmic ray muons, whose rate is 
independent of LHC luminosity, dominate the data samples 
before the muon segment veto is applied. 

The reconstruction efficiency is modified by the dif- 
ferent stopping models, since they affect the stopping 
locations of the i?-hadrons in the calorimeter. Both 
channels' signal reconstruction efficiencies are summarised 
in Table |3l Very energetic i?-hadron decays are more 
likely to be reconstructed as several jets, due to their 
large spatial extent. This causes the single jet channel 
efficiency to drop slightly, while the multi-jet channel 
efficiency grows, for heavier i?-hadrons. 



7 Systematic uncertainties on the signal 

A variety of systematic uncertainties are investigated; 
their descriptions and magnitudes are given in the fol- 
lowing list. Adding these uncertainties in quadrature 



yields a total systematic uncertainty of 23% on the pre- 
dicted number of signal events. 



— i?-hadron Nuclear Interaction (17%): To account for 
theoretical uncertainties in the interaction model, 
the cross-section for the quarks in the 7?-hadron to 
interact with a nucleus is varied up and down by a 
factor of two. The stopping fraction is recalculated 
and the largest deviation from the nominal value 
is taken as the uncertainty. The anticorrelation be- 
tween the reconstruction efficiency and the stopping 
efficiency is conservatively not taken into account. 
This procedure is applied for both stopping models, 
and the larger deviation is quoted. 

— Selection Criteria (9.9%): Each individual selection 
criterion is varied by ±10% except for the ngo re- 
quirement which is varied by zblcell. The largest 
difference across all simulated signal samples and 
selection flow is used as the systematic uncertainty. 
The 10% variation serves as a conservative estimate 
of the effects of jet energy scale and resolution un- 
certainties and other detector effects. This variation 
includes effects arising from the jets originating in- 
side the calorimeter and propagating in unusual di- 
rections. This number also accounts for the uncer- 
tainty from the limited Monte Carlo statistics. 

— Luminosity (3.4%): This uncertainty is taken di- 
rectly from Ref. [25] . 

— Calorimeter Timing (3%): In the signal simulation, 
each i?-hadron decay is given a random timing off- 
set between ±50 ns relative to the nominal bunch 
crossing time. This allows the quantification of the 
calorimeter response to significantly out-of-time de- 
cays. In reality, a signal decay will cause the trig- 
ger to be fired in the first 25 ns window in which 
it deposits enough transverse energy to satisfy the 
calorimeter trigger threshold (10 GeV); the simula- 
tion, however, only models the trigger in one bunch 
crossing. To take this into account, the signal effi- 
ciency is measured in the following windows [-10, 
15], [-5, 20] and [0, 25] ns relative to the nominal 
bunch crossing time. A 3% fractional variation in 
the efficiency is observed. The window is conserva- 
tively varied by 5 ns since the calorimeter channel- 
to-channel timing uniformity is of order 2 ns. 

— Production Cross-Section (10%): The same proce- 
dure as in Ref. [TT] is used. The production cross- 
section from PROSPINO [27j is calculated using the 
sparticle mass as the renormalisation scale with un- 
certainties estimated by varying the renormalisation 
and factorisation scales upward and downward by a 
factor of two. 
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Table 3 Signal reconstruction efficiency for various gluino mass values with a fixed x? mass of f 00 GeV. Efficiencies are shown 
for the single jet (multi-jet) channels and the uncertainties quoted are statistical only. Only the process g — > gx? i s considered 
here. These values do not take into account stopping inefficiencies, timing inefficiency, or losses from random muon segments 
discussed in Section 3] and Section 16.51 

mj (GeV) | Generic (%) \ Regge (%) 



200 


0.9±0.3 (0.1±0.1) 


1.6±0.3 (0.0±0.0) 


300 


8.8±0.6 (4.1±0.4) 


10.3±0.8 (6.1±0.6) 


400 


9.5±0.8 (6.5±0.7) 


13.1±0.9 (8.9±0.7) 


500 


8.6±0.8 (9.0±0.8) 


9.9±0.8 (12.0±0.8) 


600 


8.3±0.7 (9.8±0.8) 


9.8±0.7 (14.1±0.8) 



Table 4 The number of observed and expected events for each channel after the muon segment veto and jet energy requirement 
are applied, as in Table [2] The cosmic ray muon and beam-halo expectations are derived from the background sample and 
unpaired crossings, respectively. Only statistical uncertainties are shown in the table. 



Selection 




Number of events 




criterion 


Expected cosmic ray 


muons Expected beam-halo 


Observed 


Single Jet: Muon Veto 


4± 1 


4 ±1.8 
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Single Jet: > 100 GeV 


0.3 ±0.3 


0.5 ±0.2 





Multi-Jet: Muon Veto 


1±0.6 


0.8 ±0.8 


3 


Multi-Jet: > 100 GeV 


0.6 ±0.4 


0.2 ±0.2 


1 



8 Background estimation 

To quantify the expected number of background events 
in the search sample we investigate three different sources: 
cosmic ray muons, beam-halo muons, and calorimeter 
noise. Since cosmic ray muons occur at a constant rate 
regardless of LHC conditions their rate is measured us- 
ing the background sample. From this rate, 0.3 ± 0.3 
(0.6 ±0.4) events are expected to pass the single (multi- 
) jet selection in the search sample, as shown in Table[5J 
Here the uncertainties are statistical only. 

The rate of beam-halo muons changes with the LHC 
conditions, most importantly with the beam current. 
Beam-halo muons should contribute significantly only 
to the search sample where the beam current is far 
higher than for the other data samples. The number of 
beam-halo muons expected to pass the selection is cal- 
culated using data collected in the unpaired crossings 
from the same LHC fills as the empty crossing data. 
These events have a filled bunch in only one beam and, 
after requiring a far forward muon segment (\r]\ > 2.2), 
correspond to a nearly pure beam-halo sample. Of the 
beam-halo events that passed the single (multi-) jet cri- 
teria, all but 0.6% (0.05%) leave a muon segment. In the 
empty bunches there are 83 (531) events with far for- 
ward muon segments consistent with beam-halo that 
otherwise pass the selection. Thus 0.5±0.2 (0.2±0.2) 
beam-halo events are expected to pass the selection 
criteria without leaving a muon segment in the search 
sample. This expected yield of beam-halo background is 
added to the final numbers in the background sample of 
Table[2]to obtain the overall expected background yield 



shown in Table [H After the jet cleaning selection much 
less than 1 event is expected from calorimeter noise in 
either of the search channels. 



9 Results 



The leading jet energy distribution of the selected events 



is shown in Figures 3(a) and 3(b) for the single jet and 
multi-jet channels respectively. No excess of events is 
observed in either signal region beyond leading jet en- 
ergies of 100 GeV. 

Limits are set for the single and multi-jet measure- 
ments separately; however, since the single jet channel 
gives an expected limit on the cross-section which is 
half that of the multi-jet channel, the observed limit 
from the single jet channel is quoted as the final result. 
This limit applies to i?-hadrons decaying to a gluon and 
lightest supersymmetric particle (to^o = 100 GeV) with 
a lifetime between 10~ 5 to 10 3 seconds as described in 
Section 0] The limit takes into account systematic un- 
certainties on the signal cross-section and efficiency as 
well as on the background estimate, as described above. 

Limits are set on the signal cross-section for each 
gluino mass, mg, using a Bayesian method [28 with 
a uniform prior. Given the expected cross-section as a 
function of mass, and a limit on the expected number 
of signal events, gluino i?-hadrons with a mass mg < 
341 GeV (294 GeV) are excluded at 95% credibility level 
(C.L.) for the Generic (Regge) model. Gluino masses 
mg < 200 GeV are not investigated. Figure H] shows the 
expected and observed 95% C.L. limit and the ±1 and 
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Fig. 3 Distributions of jet energy for the single jet (a) and 
multi-jet (b) channels, for selected events before the leading 
jet energy requirement is applied. In each case the signal re- 
gion is defined as those events which have leading jet energy 
above 100 GeV. The background and control samples predict 
cosmic-ray muon yields, but do not include the background 
from beam halo. The signal distributions shown correspond 
to gluino decays to a gluon and a X°with 100 GeV mass. The 
signal distributions do not change for any gluino fl-hadron 
lifetime between 10 -5 and 10 3 seconds. 



±2 standard deviation bands. The cross-section limit 
changes rapidly in the 200 < mg < 300 GeV region as, 
with the increasing gluino mass, more decays pass the 
100 GeV jet energy selection. Above 300 GeV the limit 
depends only moderately on the gluino mass. The — 2a 
and —lcr variations of the expected limit coincide with 
the observed limit since no events are observed in the 
final selected data sample, and the background is low. 

Both the DO and CMS collaborations performed 
searches yielding comparable results for similar lifetime 
ranges. The DO collaboration analysed 410 pb~ 1 of proton- 
antiproton data collected from collisions at = 1.96 TeV 
and derived a 95% C.L. limit, excluding rrig < 270 GeV. 
The CMS collaboration used 10pb _1 of proton-proton 



200 250 300 350 400 450 500 550 600 

M a (GeV) 

Fig. 4 The 95% C.L. limit on the gluino pair production 
cross-section as a function of gluino mass, rrig, assuming 
m^,o=100GeV. Limits are shown for the Generic stopping 

model and g — > gX° decay. 



data collected from collisions at ^/s = 7 TeV. They de- 
rived a 95% C.L. limit, excluding m g < 370 GeV. 



10 Summary 

A search is presented for long-lived gluinos which have 
stopped in the ATLAS detector, using 7 TeV proton- 
proton collisions at the LHC. No evidence for the sub- 
sequent decay of these particles into gX° is found, in 
a dataset with peak instantaneous luminosity of 2 x 
10 32 cm -2 s _1 and an integrated luminosity of 31 pb -1 , 
during time periods where there are no proton-proton 
collisions. A dedicated calorimeter trigger is employed, 
and the observed events are consistent with the back- 
ground expectation derived from control samples in data. 
Limits on the gluino pair production as a function of 
gluino lifetime are derived and exclude 200 < rrig < 
341 GeV at 95% C.L. for lifetimes between 10~ 5 to 10 3 
seconds with a fixed ncutralino mass of m<>o = 100GeV 
and the Generic matter interaction model. 
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